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Abstract 

We present a new input parameter set of the Pagel model (Pagel & Tautvaisiene 1998) for the Large 
Magellanic Cloud (LMC) in order to reproduce the observations, including the star formation rate (SFR) 
history. It is concluded that the probability for (3-8)M Q stars to explode as SNe la has to be quite high 
(~ 0.17) in the LMC. As a result, a steep initial mass function (IMF) slope and existence of the outflow are 
not needed in order to attain the low [O/Fe] ratio in the LMC. As for the current supernova ratio, a high 
ratio (~ 1.3) is concluded by the new parameter set, which is consistent with the recent X-ray observations. 

Key words: galaxies: chemical evolution — Magellanic Clouds — Large Magellanic Cloud — stars: 
abundances — supernovae: general 

1. Introduction 

Chemical evolution model for a galaxy is a tool from which we can infer the history of the galaxy. A chemical 
evolution model will be believed to be right if it can reproduce the observations about the galaxy, such as the 
age-metallicity relation, the present Fe distribution of stars, and the [O/Fe] versus [Fe/H] diagram. As a result, 
we can infer some important informations about the events in the galaxy, such as the inflow rate of the material 
and the chemical compositions of Type Ia/II supernovae (SNe Ia/SNe II) using such a 'right' chemical evolution 
model. These informations give some constraints on other theories and numerical calculations. For example, people 
who are devoted to the formation of a galaxy or the nucleosynthesis in a supernova should take these constraints 
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into consideration. If they can meet them, it means that their calculations are supported by the chemical evolution 

model. This is the reason why we try to construct a good chemical evolution model. 

There are a lot of parameters in a chemical evolution model. We can not know a priori whether the number of 
the parameters is too much or too little. One concerned with a chemical evolution model should try to reduce the 
number of the free parameters in his model, while he should try to reproduce the observational data as many as 
possible. Unless such an effort, degenerate solutions that can explain some 'selected' observations will be reported 
one after another to the world. Even worse, if the history of the galaxy can not be determined uniquely using these 
solutions, we can not extract important informations about the events in the galaxy. This means that the aim of 
the chemical evolution model is not accomplished at all. 

In this study, we will consider the chemical evolution in the Large Magellanic Cloud (LMC). There are some 
excellent works about it. For example, Tsujimoto et al (1995) and Pagel and Tautvaisiene (1998) reproduce the 
age-metallicity relation, the [0/Fe]-[Fe/H] diagram, and the present mass fraction of gas very well. However, some 
of the conclusions presented in Tsujimoto et al (1995) are incompatible with those in Pagel and Tautvaisiene (1998). 
The Tsujimoto model predicts (i) a steeper initial mass function (IMF) slope (~ 1.55 - 1.72) compared to that in 
our galaxy (~ 1.35; Salpcter 1955) and (ii) a higher relative frequency of Type la to Type II supernovae (Ni^/Nu ~ 
0.2 - 0.25) compared to that in our galaxy (0.15; Tsujimoto et al. 1995). On the other hand, the Pagel model does 
not require (i) a steeper IMF slope and does not predict (ii) a higher Ni a /Nn compared to those in our galaxy. The 
concluded histories of the star formation and the inflow rate of the material are also different between their models. 
This means that their solutions are the degenerate ones mentioned above. This situation is contrast to that for our 
galaxy. Their solutions for our galaxy resemble each other (Tsujimoto et al. 1995; Pagel & Tautvaisiene 1995). 

It is true that we can not always obtain an unique solution even if we try to. Various solutions will be allowed 
when no strict constraint is derived from the observations. In fact, the number of the observations for the LMC 
is fewer than that for our galaxy. However, we think they do not make full use of the observational data for the 
LMC. For example, observations of the present Fe distribution of stars which will reflect the star formation rate 
(SFR) history are used to determine the parameters for the model of our galaxy (Tsujimoto et al. 1995; Pagel & 
Tautvaisiene 1995), while such observations are not used for the model of the LMC (Tsujimoto et al. 1995; Pagel 
& Tautvaisiene 1998). Even worse, their works contain trivial inconsistency in their analysis (see subsection 2.2 for 
details). That is why we should obtain a solution for the LMC using the observations of the SFR and a consistent 
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analysis before we conclude that the history of the LMC can not be determined uniquely by a chemical evolution 
model. 

In this paper, we use the Pagel model and determine the values of the free parameters in his model by the 
observations of the LMC, including the observation of the SFR history (Olszewski 1993). We will discuss the 
meaning of the predictions derived from the model doing a consistent analysis. In section 2, formulation of the Pagel 
model is explained. Results are shown in section 3. Summary and discussion are presented in section 4. 

2. Formulations 

2.1. The Pagel Model 

In this subsection, we present the formulation of the Pagel model. We consider one zone and assume that the gas 
is distributed uniformly and the heavy elements are well-mixed within the zone. Once we make such an assumption, 
any physical quantum is determined uniquely as a function of time in the model. This means that we do not try 
to reproduce the dispersions of the observations from the beginning. The aim of such a model is to reproduce the 
averaged value of the observations. If it can not, it will lose its justification for existence. 

They introduced a dimensionless time-like variable u, defined by 



Jo 

where w(t) is the transition probability for diffuse material ('gas') to change into stars in unit time at time t. This 
means the simple linear star formation law 



where s(u) — aS(u) is the mass in the form of long-lived stars (including compact remnants), and S is the mass 
of all stars that were ever born up to time t. a is the lock-up fraction (assumed constant). The fate of a star is 
divided into three classes; a long-lived star, Type la supernova, and Type II supernova. The ratio of their occurrence 
possibility is assumed to be a:a\\l-a-a\. g{u) is the mass of gas in the system under consideration. 
The evolution of the total baryon mass of the system is given by 




(1) 



ds ds 



(2) 



dt =UJ9 or Tu =g 



(3) 
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where f(t) and e(f) are the inflow and outflow rates of the material, respectively. Inflow is assumed to occur at a 

rate 

f(t)=u(t)e~\ (4) 

The inflowing material is assumed to be unprocessed. On the other hand, outflow is assumed to occur at a rate 
ds 

e(i) = 11— = r/ug. (5) 
T) is a free parameter. 

We explain the evolution of the metallicity. The mass fraction of ith element in the gas, Z i} is divided into two 
terms. One is the term for the instantaneous recycling, and the other is the term for the delayed production, 
%2,i- Type Il/Ia supernovae contribute to the evolution of and Zi.ii respectively. The life-time of Type II 
supernova is much shorter (~ several Myrs) than the age of the LMC (<~ 15 Gyrs). So Type II supernovae can 
be treated to explode when their progenitors are born in a chemical evolution model. This is the meaning of the 
instantaneous recycling. On the other hand, the life-time of Type la supernovae (few Gyrs) is comparable to the 
age of the LMC. So their life-time is taken into consideration in a chemical evolution model. This is the meaning of 
the delayed production. The evolution of Z\^ obeys the following equation: 

-( 9 Z hi ) =- { l-a-a 1 )x i --aZ 1 , i {*i = ^) (6) 

ds 1 — a — a\ ds 
at a at 

= ujgPis - ujgZi.i 

where m^n and Mu are averaged mass of ith heavy element and total mass of Type II supernova, m^u is defined as 

f™" m, u(m)(j)(i r n)m^ 1 dm 

™>i,n = — ' f mL , — n with = m x . (7) 

J m[ (p(m)m~ L dm 

mi^i(m) is the ith heavy-element mass produced in a star of main-sequence mass to. 4>{m) is the initial mass 
function (IMF) and x is the slope of the IMF. The bounds for SNe II progenitors were taken to be mi = 1OM and 

m u = 50M©, respectively (Tsujimoto et al. 1995). Mu is defined as 

r m " d>(m)dm 
J mi <p\m)m 'dm 

The evolution of Z 2 a is given by 

d dSjt - tu) dS ( m iM \ 

di {9Z ^ = dt am -di aZ ^ v^i^-j (9) 



(10) 
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ds(t - ii a ) a x ds 

— r, Vi - ~rr^2,i 

dt a dt 

ds 

= W(* - ti a )g(t - ii a )P 2 , 4 - -77 Z 2,i 

dt 

where t\ a is the averaged lifetime of SNe la progenitors. 

In their formulation, the frequency of SNe la ever occurred relative to SNe II is obtained as 

Nia_ _ P2,i mj,n s(u now - uti a ) 
N u Pi.i m it i a s(u 

now ) 

where u n ow is the present value for u. When we ignore the lifetime of supernova progenitors, Ni a /Nu is also written 
as 

at A f 8M ° m~( 1 + x '> 

Ni a _ J3M0 f s 

where A is the probability for (3-8)M Q stars to explode as SNe la (Tsujimoto et al. 1995). So, we can estimate the 
value of A when Ni a /Nu and x are determined from a chemical evolution model. 
As for the current supernova rates, we can calculate from Eq. (|J) and ([n]) as 

Ni a _ P 2 ,i rriiji u>(t - k a )g(t - k a ) ^ 

We also present here another physical quantum which is used in our analysis. That is r, which is the contribution 
of SNe la to the enrichment of heavy-elements in the gas. r is defined as 

r = (13) 

LUi a Mi a Ni a + LuuMiiNii 

where uii a and ton are mass fractions of heavy-element ejected into the interstellar gas from SNe la and SNe II, 
respectively. uji a and wn are defined as 

w ia = ^ : ^ : y (14) 

C s S^s,Fc + C g gZ gi Fc + Cout-^out.Fe 

u)\\ = — — — - (15) 

sZ Si o + gZ gt o + m ou t^out,o 

where Z g are the heavy-element abundance in unit mass of the gas. Z StOUt are the heavy-element abundances averaged 
over the mctallicity distribution of stars and outflow, respectively. m ou t is the total mass that is ejected from the 
system by the outflow. The factor c s , g , ut is introduced to correct the non-negligible SNe II contribution in the iron 
abundance. These are represented as 

Cg , s , out = 1 - io^ y^ , ( i6) 
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Using Eq. (|l0|), r is also represented as 

(17) 



1 , Pl.i Mil T "i,Ia Mil s("now) 

f*2 i i^Ia m i,II ^Ia s ( u now ~^^Ia) 



The i means that ith heavy element is used in estimating r. We emphasize that r does not depend on i when results 
of the calculation of supernova nucleosynthesis are used. This is because P\^ and P%^ are proportional to miji and 
m^ia, respectively (see Eq. ([j|) and (g)). 

We can also infer the value r using the present abundance pattern (Russel & Dopita 1992) and results of supernova 
nucleosynthesis as follows (Tsujimoto et al. 1995). We define the abundance pattern Xi to be compared with the 
i.,lmc as 

/ x m iJa. . /, .mill 

Xi( - r) = r l^ + il - r) M^> (18) 

and the most probable value of r — r p is determined by minimizing the following function (Yanagida et al. 1990): 

n 

g{r) = ^[loga^LMC - \ogXi{r)] 2 /n. (19) 

i=l 

Tsujimoto et al. (1995) concluded that r p — 0.16 for the LMC from this fitting (see also Figure 3) using 10 elements 
(O, Ne, Mg, S, Ar, Ca, Cr, Mn, Fe, and Ni). 

2.2. Our Standpoint on Determining the Parameters in the Paget Model 



In this subsection, we explain our standpoint on determining the parameters in the Pagcl Model. 

Pagel & Tautvaisiene (1998) determined the parameters in their model so as to explain the age-metallicity relation, 
the present mass fraction of gas, and [a/Fe] versus [Fe/H] diagrams in the LMC. a represents a-nuclei, such as O, 
Mg, Si, Ca, and Ti. 

In this study, we add two observational facts, because we should try to reproduce the observational data as many 
as possible. One is the present Fe distribution of stars in the LMC (Olszewski 1993). It will reflect the SFR history. 
The other is the contribution of SNe la to the enrichment of heavy-elements in the gas, r. 

At first, we can determine the parameters for the gas dynamics, such as g(0), s(0), co(t), tj a , and r\ as well as P\^ e 
and i"2,Fc so as to reproduce the observations of the age-metallicity relation, the present mass fraction of gas, and 
the SFR history. 
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Next, [a/Fe] versus [Fe/H] diagrams are used to extract the contribution of SNe II because a-nuclei, especially 

oxygen, are synthesized mainly in SNe II (Hashimoto 1995). This means that P l ct and P 2 ,a are determined by using 

the observations of [a/Fe] versus [Fe/H] diagrams. In this study, the [O/Fe] versus [Fc/H] diagram is mainly used to 

extract SNe IPs contribution. It is also noted that oxygen is synthesized at the outer layer in a star. So the amount 

of oxygen in the supernova ejecta is not influenced so much by the location of the mass cut (the boundary between 

the ejecta and the central compact object). On the other hand, most of the other a-nuclei mentioned above are 

synthesized near the location of the mass cut, which means their amounts in the ejecta are sensitive to the location 

of the mass cut. That is, the uncertainty of the calculations of supernova nucleosynthesis will be smallest when the 

amount of oxygen is considered. Of course, we also discuss the dependence of the results on the a-nuclei to be used 

in this analysis. This means to check the validity of the calculation of the supernova nucleosynthesis itself. 

At this point, we have to explain the problems with P\i and Pi,i- As is clear from Eq. (^|) and (j^), P\^ an d Pi,i 

are proportional to m^n and TOjjaj respectively. For example, -Pi,fo, P\,Oi ^.Fc, and Pi,o meet the relation 

EhO^H^lL and P±o = I2om^ (20) 

Pl,Fe WlFe,II P2,Fc ?™Fc,Ia 

We use the results of calculations for SNe la (Nomoto et al. 1984). So the ratio P2,o/P2,Fe is fixed from the 
beginning. In other words, we have only one parameter, Pi,o, in order to fit the [a/Fe] versus [Fe/H] diagrams. If 
we can not reproduce the observations of the [a/Fe] versus [Fe/H] diagrams completely, it reflects the problem with 
the calculation of the nucleosynthesis in a supernova explosion. 

However, Pagel and Tautvaisiene (1998) regarded all of P\^ and Pa,i to be free parameters. At the same time, 
they used Eq. (|l0|) and Eq. (|l2|) in order to estimate the occurrence frequency of SNe la relative to SNe II. It 
is quite apparent that these treatments are incompatible with each other. Moreover, they artificially shifted data 
points upwards by 0.2 dex for oxygen and downwards by 0.1 dex for silicon so as to coincide with the theoretical 
curves. This means that observational data points and theoretical curves are shifted freely and we can not extract 
useful informations from the chemical evolution model for the LMC. We emphasize again that we have to determine 
what is assumed and what is obtained in the analysis of the chemical evolution model. In this study, we believe the 
observational data points and calculation of the supernova nucleosynthesis. That is, these data points are not shifted 
artificially and results of the calculations are not changed artificially. Parameters in the chemical evolution model 
are determined using these informations. This is our standpoint on determining the parameters in the chemical 
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evolution model. 

Finally, the IMF slope is determined from Eq. (0) and A is determined from Eq. ( |Io| ) and Eq. (11) using the results 
of calculations for SNe II (Hashimoto 1995). We can also estimate the occurrence frequency of SNe la relative to 
SNe II using Eq. @ and Eq. @. 



Results 



3.1. Determination of the Parameters 



Values of the input parameters are tabulated in Table 1. Those used in the Pagel model and used in this study 
(New model) are presented in the table. P\^ and Pi,i are written in units of solar abundance of the corresponding 
element. Tq and ti a are the age of the LMC and the averaged lifetime of SNe la. uj and time are written in units of 
Gyr -1 and Gyr, respectively. g(t = Q)/m(t = To) is the ratio of the initial mass of gas to the final (t — To) total 
baryon mass. s(t — 0) is the initial mass in the form of stars. 

The comparison between the theory and the observations is presented in Figure 1 and Figure 2. In Figure 1, 
theoretical curves of the Pagel model with observations of the LMC are presented. The age-gas fraction relation, 
the age-metallicity relation, the [O/Fe] versus [Fe/H] diagram, and the present Fe distribution of long-lived stars are 
shown, respectively. In Figure 2, theoretical curves with new input parameters are shown. It is noted that Pagel and 
Tautvaisiene (1998) shifted data points for [O/Fe] upwards by 0.2 dex in their paper. In this study, such a treatment 
is not done. In the [O/Fe] versus [Fe/H] diagram, theoretical curves for Pi^o/Zq — 0.0644 (solid curve) and Pi.o/Zq 
= 0.154 (short-dashed curve) are shown for comparison. We can find clearly that solid curve reproduces better the 
[O/Fe] versus [Fe/H] diagram. So we take Pi } o/Zq — 0.0644 as a new parameter. 

Values of the output parameters are shown in Table 2. wi a , luu 7 Ni^/Nu, 7Vi a /7Vn, r, x, and A are shown in the 
table. \ 2 is the value of the chi-square probability function when data sets for the SFR are used. | [O/Fe] — [0/Fe] o b s | 2 
is the mean square of the difference between the observed [O/Fe] and calculated one at the same [Fe/H]. Data points 
for the [O/Fe] are shifted by 0.2 dex when the value of the |[0/Fc] — [0/Fe] o b s | 2 for the Pagel model is calculated. 
As for r, the minimizing function g(r) is shown in Figure 3. r = 0.16 (solid vertical line) is the most probable value 
from the fitting of Eq. (|l9|). r = 0.12 and 0.20 are the calculated values by the Pagal and New models, respectively. 
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Finally, time evolution of the total mass (solid curve), mass of gas (short-dashed curve), and mass of stars (long- 
dashed curve) concluded by the Pagel and New models arc presented in Figure 4. Time is defined as 14 - Age 
[Gyr]. 

3.2. Problem with the Supernova Nucleosynthesis 

In this subsection, we compare the theoretical curves for the [a/Fe] versus [Fe/H] diagrams with the observational 
data points, a represents a-nuclei, such as Mg, Si, Ca, and Ti. We can check the validity of the calculation of the 
supernova nucleosynthesis due to this test. We emphasize again that we have no free parameters already, because 
fi,i=M g ,Si,Ca,Ti are determined uniquely when Pi,Fe and P^o are given. 

The results are shown in Figure 5. We can find that [Mg/Fe], [Si/Fe], and [Ca/Fe] are not well reproduced as 
long as Px^o/Zq = 0.0644 is adopted. On the other hand, [O/Fe] is not well reproduced when Px^o/Zq = 0.154 is 
adopted (see Figure 2). This is the problem what Pagel & Tautvaisiene (1998) pointed out. 

As for the ratio of [Ti/Fe], both of the theoretical curves can not reproduce the data points well. That is, more 
of Ti is required in order to explain the observations. This shortage problem of Ti will mean the shortage problem 
with the supernova nucleosynthesis (Tsujimoto et al. 1995; Hashimoto 1995). 

Finally, the output parameters for Pi^/Zq = 0.154 are shown in Table 3. The uncertainties of the output 
parameters will be inferred when we compare the values in Table 3 with those in Table 2. 

4. Summary and Discussion 

As stated in section 1, chemical evolution model for a galaxy is a tool from which we can infer the history of the 
galaxy. We also emphasize again that a chemical evolution model has to contain many parameters in it. We have to 
determine these parameters using results of calculations of supernova nucleosynthesis and observations in a galaxy. 
When these results and/or observations are shifted artificially, as done in Pagel & Tautvaisiene (1998), almost all 
informations arc treated as free parameters. Wc will be able to extract very little informations from such works. We 
must always keep in mind what is assumed and what is derived in the study of the chemical evolution in a galaxy In 
this study, observations and results of calculations of supernova nucleosynthesis are believed and used to determine 
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the parameters in the chemical evolution model. In this section, we discuss what the New model presented in this 
study suggests as a result. 

As for the ratio g(t — 0)/m(t — Tq), the New model requires relatively high ratio (see Table 1 and Figure 4). 
This means that we require an initial condition in which star formation is forbidden and only gas exists. What 
can suppress the star formation at zero metallicity in the LMC? In a simplest situation, a cloud with only thermal 
support, collapse should occur if the mass exceeds the Jeans (1928) mass, 



where Tk is the kinetic temperature (K), p is the mass density (g cm -3 ), and n is the total particle density (cm -3 ). 
So, low temperature and/or high number density are required in order to form a star or, at least, to form a molecular 
gas cloud. So, it is suggested that the LMC was born in which the temperature is high and/or the particle density 
is low enough to prevent gas from collapsing and from forming a molecular gas cloud. 

We discuss the effect of the outflow. Pagel & Tautvaisiene (1998) introduced an outflow in order to reduce the 
metallicity of the system. The [O/Fe] ratio is also reduced by the effect of the outflow when a star burst phase is 
assumed. This is because the material of SNe II is ejected more than that of SNe la. However, the Pagel model can 
not reproduce the SFR history very well (see Figure 1 and the value of the x 2 probability function in Table 2). On 
the other hand, the New model can reproduce the present low metallicity and the SFR history at the same time 
although the New model requires very smaller value for rj than the Pagel model. Moreover, the [O/Fe] versus [Fe/H] 
diagram can be reproduced well without shifting the data points since a low [O/Fe] ratio can be attained by using 
a high value of A, not by using the effect of the outflow. We could not reproduce the observations well when r\ is 
set to be 1.0 like the Pagel model. Although we can not conclude that our solution is an unique one, one has to 
construct a new parameter set which can reproduce the observations as well as ours if he insists the existence of the 
outflow in the LMC. Additionally, a steeper IMF slope (Tsujimoto et al. 1995) is not needed in order to attain the 
low value of the [O/Fe] since a high value of A can realize it. 

The ratio Ni a /Nn inferred from the New model is very high (=0.48) compared with other models (Tsujimoto et 
al. 1995; see also Table 2). This is because high ratio of P2,Fe/ Pi,Fe (i-e. high ratio of A) is assumed in the New 
model (see Table 1). However, r=0.20 in the New model is close to that inferred from the observations (r=0.16; see 
Eq. (|l9|)) when the weak dependence of g(r) on r around r ~ 0.16 is considered (see Figure 3). That is why we 




(21) 
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think the high ratio of Ni a /Nu is not ruled out from the present chemical composition in the LMC. As for the ratio 

Nia/Nu, a higher ratio (= 1.3) is concluded by the New model than the Pagel model. Since such a high ratio, of 

order 1, is estimated by X-ray observations (Hughes et al. 1995), the New model is thought to be consistent with 

the observations. 

We found that the problems with the supernova nucleosynthesis using the [a/Fe] versus [Fe/H] diagrams. We can 
not solve these problems now. However, these problems may be solved if the metallicity effect for the supernova 
nucleosynthesis are taken into consideration (Woosley & Weaver 1995). We will perform such calculations and report 
its result in the near future. As for the shortage problem of Ti (Tsujimoto et al. 1995; Hashimoto 1995), asymmetric 
explosion models for SNe II may solve the problem (Nagataki 1997, 1999). This is because more of Ti is produced in 
the asymmetric explosion model by the strong alpha-rich freezeout (Nagataki 1997, 1999). We will discuss whether 
these models can explain the [a/Fe] versus [Fe/H] diagrams at the same time in the forth coming paper. 

We present a new input parameter set of the Pagel model which reproduces the observations of the LMC very 
well. The aim of constructing a chemical evolution model is to give some constraints on other theories and numerical 
calculations. In order to attain such an aim, the chemical evolution model has to explain the observations in the 
galaxy as many as possible. At the same time, must always keep in mind what is assumed and what is derived in 
the study of the chemical evolution in a galaxy. Since the parameters in a chemical evolution model are determined 
by the observations and results of calculations of supernova nucleosynthesis, we will be able to determine these 
parameters more precisely by further precise observations and/or calculations. We hope more observations and 
further discussions will be presented by many groups in the world in order to understand the system of the Large 
Magellanic Cloud. 

We are grateful to Pagel, B.E.J, for useful comments. This research has been supported in part by a Grant-in- Aid 
for the Center-of-Excellence (COE) Research (07CE2002) and for the Scientific Research Fund (7449, 199908802) of 
the Ministry of Education, Science, Sports and Culture in Japan and by Japan Society for the Promotion of Science 
Postdoctoral Fellowships for Research Abroad. 
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Table 1. Values of the Input Parameters. 





Pagel Model 


New Model 




0.28 


0.07 


Pi,o/Z° 


0.70 


0.0644 


-F^Fc/Zq 6 


0.42 


0.25 


T g (Gyr) 


14 


14 


iia (Gyr) 


1.33 


2.2 


V 


1.0 


0.001 


Lu(t) 


0.15 (t < 2) 


0.35 {t < 0.2) 




0.08 (2 < t < 11) 


0.12 (0.2 < t < 11) 




0.50 (11 < t < 14) 


0.45 (11 < t < 14) 


g(0)/m(T G ) 





0.82 


s(t = 0) 









Left column: those used in the Pagel model. Right: those used in this study. and P 2 ,i are written in units of 
solar abundance of the corresponding element. Tq and t\ a are the age of the LMC and the averaged lifetime of SNe 
la. oj and time are written in units of Gyr -1 and Gyr, respectively. g(t = 0)/m(t = Tq) is the ratio of the initial 
mass of gas to the final (t = Tq) total baryon mass. s(t = 0) is the initial mass in the form of stars. 
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Table 2. Values of the Output Parameters. 





-TagCl ivlOQCl 


INcW IVlOUCi 


Wla 


(J. 21 


not 

0.25 




0.14 


0.21 


N Ia /N u 


0.20 


0.48 


Nu/Nn 


0.42 


1.3 


r 


0.12 


0.20 


X 


1.01 


1.06 


A 


0.075 


0.17 


x 2 


0.41 


0.73 


|[0/Fe]-[0/Fe] obs | 2 


0.19 


0.19 



Left column: those derived by the Pagel model. Right: those derived by this study. u>\ a and u>n are mass fractions 
of heavy-element ejected into the interstellar gas from SNe la and SNe II, respectively. N^/Nu and Ni a /Nu are the 
occurrence frequency of SNe la relative to SNe II and the ratio of the current supernova rates, respectively, r is the 
contribution of SNe la to the enrichment of heavy-elements in the gas. x is the IMF slope. A is the probability for 
3-8M stars to explode as SNe la. x 2 is the value of the chi-square probability function when data sets for the star 
formation rate are used. |[0/Fe] — [0/Fe] o bs| 2 is the mean square of the difference between the observed [O/Fe] and 
calculated one at the same [Fe/H]. 
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Table 3. Values of the Output Parameters for Pi,o/Zg = 0.154. 





0.25 




0.21 




0.43 


Nu/Nn 


1.16 


r 


0.13 


X 


0.22 


A 


0.51 


x 2 


0.73 


|[0/Fe]-[0/Fe] obs | 2 


0.40 



Same as Table 2 but for Pi,o/Z§ = 0.154 model. 



Figure Captions 

Fig. 1. Theoretical curves of the Pagel model with observations of the LMC. The upper panel on the left: the 
age-gas fraction relation. Data point is from Westerlund (1997). The lower panel on the left: the age-metallicity 
relation. Data sources are as follows: open triangles, Olszewski et al. (1991); filled triangles, Girardi et al. 
(1995); open hexagons, Geisler et al. (1997); open squares, de Freitas Pacheco et al. (1998). The upper panel 
on the right: the [O/Fe] versus [Fe/H] diagram. Data sources are as follows: open triangles, Luck & Lambert 
(1992); open squares, Thevenin (1997); open hexagons, Hill, Andriewski & Spite (1995); filled triangles, Spite, 
Barbuy & Spite (1993); filled hexagons, Jutter et al. (1992); five — pointed stars, McWilliam & Williams 
(1991); three — pointed stars, Barbuy, de Freitas Pacheco & Castro (1994); cross, Richtler, Spite & Spite (1989); 
asterisks, Nissen & Schuster (1997). The lower panel on the right: the present Fe distribution of long-lived stars. 
Data source is from Olszewski (1993). 

Fig. 2. Same as Fig.l but with new input parameters. In the [O/Fe] versus [Fe/H] diagram, theoretical curves for 
P\,o/Zq — 0.0644 (solid curve) and Pi^q/Zq = 0.154 (short-dashed curve) are shown for comparison. 

Fig. 3. The minimizing function g(r) as a function of r. 10 elements (O, Ne, Mg, S, Ar, Ca, Cr, Mn, Fe, and Ni) are 
used for the fitting, r = 0.16 (solid vertical line) is the most probable value, r = 0.12 and 0.20 are the calculated 
values by the Pagal and New models, respectively. 

Fig. 4. Time evolution of the total mass (solid curve), mass of gas (short-dashed curve), and mass of stars (long- 
dashed curve). Time is defined as 14 - Age [Gyr]. 

Fig. 5. [a/Fe] versus [Fe/H] diagrams with the observational data points. Solid curves represent the theoretical 
curves for Piq/Zq = 0.0644. Short-dashed curves represent the ones for Pi^q/Zq = 0.154. 
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Fig. 2. 
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Fig. 4. 
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